PASJ: Publ. Astron. Soc. Japan , 1-??, 
© 2012. Astronomical Society of Japan. 



Detection of Change in Supercycles in ER Ursae Majoris 
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Abstract 

We examined data on about 20 years of observations of ER UMa available in AAVSO, VSNET, AFOEV, 
NSVS, VSOLJ databases together with published light curves. The obtained O — C diagram revealed a 
systematic change of the supercycle (time interval between two successive superotbursts) within 43.6 and 
59.2 d. The time-scale of this change variation is from 300 to ~ 1900 d. The number of normal outbursts 
within the supercycles also varied from 4 to 6 although no strong correlation between this number and 
supercycle length was found. We suggest that appearance of negative superhumps is responsible for the 
observed variations of number of normal outbursts. Our results generally confirms the expectations by the 
thermal-tidal instability theory. 

Key words: accretion, accretion disks — stars: dwarf novae — stars: individual (ER UMa) — stars: 
novae, cataclysmic variables 



1. Introduction 

Dwarf novae (DNe) are a class of cataclysmic variables 
(CVs) — close binary systems consisting of a white dwarf 
as a primary component and a red or brown-dwarf sec- 
ondary transferring matter via the Roche-lobe overflow. 
Because of the instabilities that take place in the accre- 
tion disks normal outbursts and less frequent superout- 
bursts with larger amplitudes take place (see Warner 1995, 
Hellier 2001). DNe that show both normal outbursts and 
superutbursts are called the SU UMa-type stars. These 
stars have one distinguishing property — light-curve mod- 
ulation with the periods close to, but different from, the 
orbital one that appear near the superoutburst maxima 
and sometimes persist at other times (Vogt 1980). These 
modulations are called positive superhumps if their pe- 
riod is longer then the orbital one and negative if shorter. 
Positive superhumps are believed to arise from the peri- 
odic viscous dissipation in a disk driven to resonant os- 
cillation by the tidal field of the secondary star, and pro- 
cessing slowly in the prograde direction (see Whitehurst 
1988, Hirose, Osaki 1990, Lubow 1991, Smith et al. 2007, 
Wood et al. 2011) and negative superhumps arise from 
the precession of tilted accretion disk, respectively (Wood 
et al. 2009, Montgomery, Martin 2010). 

SU UMa stars have two extreme subgroups with very 
long (about decades) and short (tens of days) supercycles. 
They are WZ Sge and ER UMa subtypes, respectively. 
There are five ER UMa-type stars known by now: RZ 
LMi (Nogami et al. 1995), DI UMa (Kato et al. 1996), 
ER UMa (Kato, Kunjaya 1995), VI 159 Ori (Robertson 
et al. 1995), IX Dra (Ishioka et al. 2001) and all of them 
show extremely short supercycle (18.9 - 54 d) and posi- 



tive superhumps (Olech et al. 2008). Moreover, ER UMa- 
type stars were known to show a pronounced stability of 
the supercycle length and outburst pattern, which can be 
particularly seen in the folded light-curves and O — C dia- 
grams presented in (Robertson et al. 1995). This stability 
and short supercycles can be explained as a result of a very 
high and constant mass-transfer rate in the framework of 
thermal tidal instability (TTI) model (Osaki 1995a). 

However, long-term monitoring of V1159 Ori revealed 
secular changes in the supercycle and a hint of the peri- 
odic behavior of such variations (Kato 2001). The TTI 
model suggests that this phenomenon is the consequence 
of a change in mass-transfer rate in the binary system, 
although the origin of these changes is still an open ques- 
tion. The hibernation scenario, or the theory of cyclic 
evolution, suggested by Livio (1992), also assumes that 
mass-transfer rate in CVs can vary secularly and proposes 
irradiation of the secondary by the white dwarf which 
is still hot after the nova eruption as one of the possi- 
ble mechanisms of an increased mass-transfer. However, 
the time-scale of M variation in this scenario is very long 
comparing to these observations. On the other hand for 
V1159 Ori, Kato (2001) proposed that solar-type cycles of 
the secondary which are possibly observed in cataclysmic 
variables (Bianchini 1988, Ak et al. 2001) may be respon- 
sible for the phenomenon. 

Here we report that another one ER UMa-type star — 
ER UMa itself, also shows variations of the supercycle. A 
discussion on the origin of this phenomenon is even more 
intriguing in the context of recent discovery by Ohshima 
et al. (2012) of persistent negative superhumps in this 
dwarf nova. The presence of negative superhumps was 
not taken into account in past studies of SU UMa-type 
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stars, while this factor may strongly affect the outburst 
behavior. Osaki, Kato (2012) have also found that the 
presence of negative superhumps is correlated with normal 
outbursts, thus it is reasonable to investigate the possibil- 
ity whether the presence of negative superhumps is related 
to the long-term variations of the outburst behavior. 

A comprehensive study of long-term variations in ER 
UMa-type stars is therefore essential for better under- 
standing the evolution of dwarf novae in general, in partic- 
ular changes in mass transfer-rate, the mechanisms that 
lead to the appearance of negative superhumps and their 
effect on the outburst behavior. 

2. Observation and Data Analysis 

2.1. Light- curves 

We analyzed all the data on ER UMa available 
in American Association of Variable Stars Observers 
(AAVSO) 1 , Variable Stars NETwork (VSNET) 2 , 
Association Frangaise des Observateurs d'Etoiles 
Variables (AFOEV) 3 , Northern Sky Variability burvey 
(NSVS) 4 and Variable Star Observers League in Japan 
(VSOLJ) 5 databases which represent more then 20 years 
of observations. The earliest superoutburst was on 1992 
April 27 (JD2448740) and the latest one on 2012 June 
12 (JD2456090.7). We also used observations obtained 
by one of the authors (S. Shugarov) in the Stara Lesna 
Observatory of Slovak Academy of Sciences, and by 
the group of E. Pavlenko in the Crimean Astrophysical 
Observatory. A long-term light-curve of ER UMa from 
these observations is presented in figure 1. 

2.2. Moments of the superoutbursts and O — C diagram 

In order to follow the supercycle change in ER UMa 
we plotted O — C diagram using the data mentioned. The 
moments of the superoutbursts and normal outbursts were 
determined visually from the light curves. The mean pe- 
riod essential for plotting such diagram was determined 
by summing up the cycle counts and was chosen to be 
48.7 d by the least-squares method. There were sev- 
eral gaps in observations usually 3-4 cycles long, partic- 
ularly around the solar conjunction, which can be also 
seen in figure 1. The time intervals between two subse- 
quent detected superoutbursts in some such cases were 
large enough to bring considerable uncertainties in calcu- 
lations of the number of cycles in such gaps. For every 
uncertain moment of the superoutburst we checked the 
two closest possible values of the cycle number, for each 
of them calculated the mean period for the few closest 
supercycles before, across and after the gap and regarded 
the variant as less probable if it leaded to the significant 
period jump. The list of the moments of superoutbursts, 
corresponding supercycle count — E (the number of su- 
peroutbursts passed since JD2448740), supercycle length, 

1 http://www.aavso.org/ 

2 http:/ /www. kusastro.kyoto-u.ac.jp/vsnet/ 

3 http:/ /cdsarc. u-strasbg.fr/afoev/ 

4 http:/ /skydot. lanl.gov/nsvs/nsvs.php 

5 http://vsolj.cetus-net.org/ 



measured by the linear fitting of successive five epochs 
of O — C diagram and the number of normal outbursts 
observed within the supercycle is presented in the table 1. 

2.3. Normal outbursts phases 

We noticed that ER UMa shows two types of outburst 
behavior with 3-4 normal outbursts within the supercycle 
observed and with 5 or more. Examples of these types 
can be seen in figure 2 and figure 3. Following Osaki, 
Kato (2012) and Smak (1985) here we use Type S for the 
stages with many normal outbursts and short intervals 
and Type L with reduced numbers of normal outbursts 
and long intervals. Because of the gaps in observations 
and because the data were not very precise for the most 
part of supercycles it was possible to determine only ap- 
proximate numbers of normal outbursts within them, or 
even only lower limits. In the table 1 we give the proba- 
ble numbers, as "3-4" for the cases were it was impossible 
to choose the unique value, and lower limits, as ">4", for 
the supercycles, where for example four normal outbursts 
were detected but the significant part of the supercycle re- 
mained unobserved, and we could not say whether there 
were other normal outbursts or not. 

3. Discussion 

3.1. O — C Behavior, Comparison with Other SU UMa 
DNe 

The O — C diagram for the whole analyzed period is pre- 
sented in figure 4. We can see that ER UMa underwent 
several stages with a relatively stable supercycle length 
(linear portion in the diagram) which lasted for 6 - 30 
cycles and terminated by an abrupt period change. The 
local mean supercycle length was determined by the linear 
fit of each time interval with relatively stable supercycle 
in the O — C diagram. It was shown that the supercycle 
varied from 43.7 at the very beginning of the O — C dia- 
gram to 59.2 d. The local values of the supercycle length 
are also presented in the table 1. 

The O — C behavior presented in figure 4 is reminiscent 
to that of V1159 Ori, shown in figure 2. in (Kato 2001). 
Although one can see such a pattern, for example between 
E = 60 and E = 125 cycles in figure 4, the time-scale of the 
variability is significantly longer in ER UMa. A detailed 
analysis of the supercycle variations in SU UMa-type stars 
was presented by (Vogt 1980). Although only ordinary SU 
UMa-stars were discussed in Vogt (1980) (we do not con- 
sider a comparison with SS Cyg-type DNe since it is out 
of the scope of this paper) , the published O — C diagrams 
have the same features as ER UMa, particularly of SU 
UMa itself (figure 1 in Vogt 1980). Supercycles of SU 
UMa remained stable for relatively long periods of time 
(10 - 20 cycles) and switched randomly between several 
characteristic values. Nevertheless, while Vogt (1980) has 
shown that the supercycles in SU UMa possess 3 charac- 
teristic values, in ER UMa there are no preferable super- 
cycle lengths as can be seen in figure 5. 
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Fig. 1. Long-term light-curve of ER UMa from AAVSO, VSNET, AFOEV, NSVS and VSOLJ observations. 
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Fig. 4. O — C diagram for ER UMa with the phases of normal outbursts and the types of superhumps observed. Open and filled 
circles correspond to L and S normal outburst phases, respectively, supercycles with undetermined phases are marked with crosses. 
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Fig. 5. Distribution of the supercycle lengths of ER UMa. 



Detection of Change in Supercycles in ER UMa 



5 




2455640 2455660 2455680 2455700 2455720 
JD 



•i • • • 

•M 

• »«*•* • • 
• • M 
• • • 



I I '"v:V . Ill 



• • • 



• • § 



• • • 
• •• # 



• ■ • 



1 1 1 

2449760 2449780 2449800 

JD 



2449820 



2449840 



Fig. 2. Example of the L type normal outburst behav- 
ior of ER UMa (March - June 2011) 

Fig. 3. Example of the S type normal outburst behav- 
ior of ER UMa (January - May 1995) 



3.2. Supercycle Length Change 

We measured the local lengths of the supercycle by the 
linear fitting of O — C diagram and plotted them against 
time. Results can be seen in figure 6 where the supercycle 
length change is presented. Supercycle lengths (T s ) ap- 
peared to change discontinuously and show oscillations 
superimposed on a secular increasing trend. The rate 
of supercycle growth is dT s /dE = 0.033(6) d/cycle and 
the period derivative is dT s /T s = 6.7(6) x 10~ 4 , respec- 
tively. We used a half-delete jackknife method described 
in (Shao, Tu 1995) to determine reliability of the obtained 
trend. It yielded that the period increase is statistically 
significant with a 90% confidence. Nevertheless, these for- 
mal error estimates should not be taken seriously. In 
Otulakowska-Hypka et al. (2012) it was shown that the 
period derivative for IX Dra is 1.8 x 10~ 3 , comparable 
with that found for ER UMa. We suppose that it is a 
typical value of supercycle change rate for ER UMa-type 
stars. In figure 5 there is a hint of bimodal distribution 
of the periods of supercycle with maxima at 48 and 53 d, 
respectively. Osaki (1995a) proposed a model of ER UMa 
type stars based on the TTI theory which successfully 
explained the light-curves and short recurrence times of 
the superoutbursts assuming extremely high mass-transfer 
rate (about 10 times higher than that expected from CV 
evolutionary scenario based on gravitational-wave radia- 
tion). It also provided a dependence between the period 
of supercycle and mass-transfer rate (presented in figure 
1 in Osaki 1995a) and the minimum possible supercycle 
length for ER UMa-type stars (about 40 d). According 
to this dependence there are two possible values of M for 
each T s . These two possibilities can be distinguished by 
considering the duty cycle of superoutbursts (ratio of the 
duration of superout burst to the supercycle length). In 
our case, the possibility of larger M can be excluded be- 
cause the observed duty cycles were much smaller than 
0.5. Observed variations in the supercycle length of ER 
UMa correspond to the change of mass-transfer rate from 
3.8 Mi6 to 2 Mi g where Mi6 is the mass-transfer rate in 
units of 10 16 g s" 1 and the increasing trend in T s means 
that the mass-transfer rate in this system is decreasing 
secularly. It is also remarkable that the most significant 
change of the supercycle and the mass-transfer rate (al- 
most twice) occurred in first break in the O — C diagram 
(JD2449932) right after the initial long stable interval. All 
the further supercycle changes were less dramatic. The pe- 
riod varied between 46.2 and 54.1 d which corresponds to 
a ~ 15 % change in mass-transfer rate. It is also notable 
that the minimum value of the supercycle was not below 
the shortest achievable supercycle in (Osaki 1995a) and we 
therefore do not need to assume an additional mechanism 
as in RZ LMi (Osaki 1995b) 

3.3. Normal Outbursts 



Another important problem in the behavior of ER UMa 
is the the correlation between the length of the supercycle 
and the number of normal outbursts within it. For the sys- 
tems with a high mass-transfer rate (like ER UMa), both 
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Fig. 6. Variation of the supercycle length measured by the linear fitting of the O — C diagram. The solid line shows the secular 
increasing trend, whose statistical significance was verified by a half-delete jackknife method (Shao, Tu 1995) with a 90% confidence. 
Open and filled circles correspond to L and S normal outburst phases, respectively, supercycles with undetermined phes marked 
with crosses. 



the waiting time for normal outbursts (T n ) and the wait- 
ing time for superoutbursts (T s ) are approximately pro- 
portional to M _1 (Osaki 1995b). In this case T n /T s does 
not depend on M and the number of normal outbursts be- 
tween two successive superoutbursts should be constant 
according to the model of ER UMa stars proposed by 
Osaki (1995a). However, as was mentioned in Subsection 
2.3, ER UMa shows two types of normal outburst phases: 
a phase with 3-4 normal outbursts observed and with more 
than 4. In a Type L phase presented in the figure 2 ER 
UMa underwent four normal outbursts (marked with ar- 
rows) that occurred once in 7 days. In figure 3 five normal 
outbursts can be seen, moreover in the second supercycle 
the 6-th normal outburst near JD2449824 was likely un- 
observed because of the gap in the observations. The time 
intervals between two subsequent normal outbursts were 
4 days. In figure 4 and figure 6, filled circles correspond to 
the supercycles with Type S normal outburst phase and 
open circles to the supercycles with Type L phase. If it 
was impossible to determine the type of the phase within 
the supercycle we mark them with crosses. It can be seen 
from these plots that there are no strong correlation be- 
tween the normal outburst phases and the length of the 
supercycle. Types L and S were observed in any value of 
the supercycle length. Bends on the O — C diagram, or 
shifts of the supercycle period, did not affect the number 
of normal outbursts. Apparently there should be some ad- 
ditional factor, which is responsible for the change of the 



normal outbursts phases. Osaki, Kato (2012), Ohshima 
et al. (2012) suggested that the type of phase may depend 
on the absence or presence of negative superhumps. This 
possibility is discussed in the following subsection. 

3.4- Negative Superhumps and Outburst Behavior 

Negative superhumps are believed to arise from the pre- 
cession of a tilted accretion disk, nevertheless the origin 
of the disk tilt is still under debate. In such cases the 
hot spot appears not only at the disks rim but deeper in 
the primaries potential, too. Montgomery, Martin (2010) 
proposed a mechanism based on the hydrodynamic lift. 
Differences between the supersonic velocities of the accre- 
tion flow under and below the disk cause the appearance 
of the lift force and make the disk tilted when the mass 
transfer is high. The accretion stream can thus reach the 
inner parts of the accretion disk. Interaction of the mag- 
netic fields with rotating plasma was also suggested as a 
possible source of the disk tilt (Murray, Armitage 1998, 
Murray et al. 2002, Smak 2009). A disk tilt is likely to af- 
fect the outburst behavior of SU UMa stars, and it is very 
important to study the correlation between the the pres- 
ence of negative superhumps (as an indicator of the tilt 
of the disks) and outbursting activity. In 2011 it became 
evident that persistent negative superhumps in ER UMa 
were observed even during the superoutburst (Ohshima 
et al. 2012) while a retrospective search for the litera- 
ture and the available data suggests that they were ab- 
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sent in the past. Kato, Kunjaya (1995) in their discovery 
paper reported the detection of positive superhumps in 
ER UMa. No negative superhumps which are considered 
as characteristic feature of ER UMa in the present state, 
were observed by them in 1994. The first possible detec- 
tion of negative superhumps in ER UMa was claimed by 
Gao et al. (1999), but the amplitudes were smaller than 
in Ohshima et al. (2012) and it is likely that, even if the 
negative superhumps did exist in 1998, they were not so 
strong in 1998 as in 2011. Since then the presence of 
the superhumps with period less than the orbital one was 
confirmed by Zhao et al. (2006), Kjurkchieva, Marchev 
(2010). In our retrospective study we also confirmed the 
presence of negative superhumps in 2010 in the AAVSO 
data. It appears likely that the negative superhumps were 
present since 2006. 

3.5. Correlation with Normal Outburst Phase 

TTI theory suggests that a mass-transfer to a tilted 
disk is expected to reduce the number of normal outbursts 
(Osaki, Kato 2012). Since the accretion stream in such 
cases can reach the inner parts of the disk, material is 
accumulating not primarily on the disks rim. It leads to 
the reduction of "outside-in" outbursts. Here we examine 
the correlation between the normal outburst phase and 
the presence of negative superhumps. In the figure 4 we 
mark the types of superhumps observed by the authors 
mentioned above. One can see that the presence of nega- 
tive superhumps seems to depend strongly on the phase of 
the normal outbursts. The first possible detection of neg- 
ative superhumps was near the E = 44 on figure 4 when 
the first L stage was observed. The next one, claimed by 
Zhao et al. (2006) was on 2004 February 29 which corre- 
sponds to the E = 88, just at the transition from S to L 
phases. Negative superhumps observed by Kjurkchieva, 
Marchev (2010) in 2010 were in E = 119 at the diagram, 
also at the beginning of L phase. According to the de- 
tailed light curves of ER UMa provided by the observa- 
tional campaign of the VSNET collaboration (Kato et al. 
2004) negative superhumps have been observed since the 
2011 (E = 135) and four normal outbursts occurred within 
each supercycle. It confirms that the phase of the normal 
outbursts can be strongly affected by the presence of neg- 
ative superhumps, and negative superhumps can indeed 
reduce the number of normal outbursts. So we can now 
assume that Type L phase can be regarded as "negative 
super hump" phase. 

3.6. Correlation with the Supercycle length 

According to figure 4, positive and negative super- 
humps were observed both on the descending and ascend- 
ing branches of the O — C diagram, and they do not cor- 
relate with the supercycle length. Negative superhumps 
were with no doubt observed since 2011. During this pe- 
riod O — C underwent a bend near E = 150 but it had 
no affect on the presence of negative superhumps. In the 
TTI model, the next superoutburst occurs when the disk 
stores a certain amount of angular momentum (J cr it), and 
supercycle is not expected to change if J cr u from the sec- 
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ondary is constant regardless of the point of accretion on 
the disk. A disk tilt therefore is not expected to affect the 
supercycle length and our analysis confirmed this expec- 
tation. 

4. Conclusions 

In this paper we examined data of about 20 years of 
observations of ER UMa. We followed the evolution of 
the change in supercycle length, the behavior of nor- 
mal outbursts and the types of superhumps observed. 
We found the oscillations of the supercycle length su- 
perimposed with a secular increasing trend and associ- 
ated it with a gradual change of the mass-transfer rate in 
the binary system. The estimated period derivative was 
dT s /T s = 6.7(6) x 10~ 4 . We also found that ER UMa shows 
two types of normal outburst phases with 3-4 and with 5 
or more normal outbursts within the supercycle. We also 
examined the possible correlation between the supercycle 
length, normal outburst phase and the presence of nega- 
tive superhumps and found a strong dependence between 
the last two. A reduced number of normal outbursts was 
observed when negative superhumps were detected. It 
confirms the assumption of Osaki, Kato (2012) that neg- 
ative superhumps tend to suppress the normal outbursts. 
No correlation between the length of the supercycle and 
other characteristics was found. 
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Table 1. Moments of the superoutbursts, supercycle lengths (T s ) and the number of normal outbursts within the supercycle (JVno)- 
Since the moments of superoutbursts were determined from several databases, multiple entries for the same E sometimes occur. 



E 


Date 


JD 


T s (d) 


iVNO 


E 


Date 


JD 


T s (d) 


N NO 





1992/04/27 


2448740 


49.6 




90 


2003/12/18 


2452992 


45.4 


3-4 


3 


1992/09/05 


2448871 


49.0 




91 


2004/01/28 


2453033 


44.0 


4 


6 


1993/01/13 


2449001 


44.1 




91 


2004/01/28 


2453033 


43.3 




7 


1993/02/27 


2449046 


44.5 




92 


2004/03/13 


2453078 


46.5 


4 


8 


1993/04/17 


2449095 


45.0 


>4 


92 


2004/03/13 


2453078 


46.8 




9 


1993/05/29 


2449137 


43.6 


>4 


93 


2004/05/01 


2453127 


47.5 


4 


10 


1993/07/12 


2449181 


AAA 




94 


2004/06/16 


2453173 


46.3 


3-4 


11 


1993/08/22 


2449222 


43.9 




94 


2004/06/16 


2453173 


46.6 




12 


1993/10/14 


2449275 


44.0 




96 


2004/09/15 


2453264 


47.2 




13 


1993/11/18 


2449310 


43.6 




97 


2004/11/05 


2453315 


47.8 




14 


1994/01/04 


2449357 


42.1 




97 


2004/11/06 


2453315 


49.2 




15 


1994/02/15 


2449399 


43.1 




98 


2004/12/25 


2453365 


48.3 


3-4 


1G 


1994/03/29 


2449441 


42.6 


>4 


99 


2005/02/10 


2453412 


49.3 


3-4 


17 


1994/05/11 


2449484 


42.5 


>4 


99 


2005/02/10 


2453412 


49.8 




18 


1994/06/24 


2449528 


42.4 


5-6 


100 


2005/04/04 


2453465 


52.1 


4 


19 


1994/08/03 


2449568 


42.1 




101 


2005/05/22 


2453513 


52.2 


3-4 


20 


1994/09/15 


2449611 


42.1 




103 


2005/09/08 


2453621 


52.4 




21 


1994/10/27 


2449653 


43.5 




104 


2005/10/29 


2453673 


53.5 




22 


1994/12/09 


2449696 


45.0 




104 


2005/10/30 


2453673 


53.4 




23 


1995/01/25 


2449743 


46.9 


>3 


105 


2005/12/23 


2453728 


54.9 




24 


1995/03/14 


2449791 


47.3 


5 


105 


2005/12/23 


2453728 


55.2 




25 


1995/05/02 


2449840 


48.71 


5 


105 


2005/12/24 


2453728 


56.0 




27 


1995/08/02 


2449932 


50.5 




106 


2006/02/17 


2453784 


53.3 




28 


1995/09/29 


2449990 


54.4 


5 


106 


2006/02/17 


2453784 


54.2 


4 


29 


1995/11/24 


2450046 


59.6 




107 


2006/04/08 


2453834 


54.7 


6 


30 


1996/01/31 


2450114 


59.4 




108 


2006/06/06 


2453893 


55.5 




31 


1996/03/25 


2450168 


55.8 


>4 


108 


2006/06/07 


2453894 


56.1 


>5 


32 


1996/05/22 


2450226 


48.9 




111 


2006/11/20 


2454060 


55.2 




33 


1996/07/04 


2450269 


47.2 




111 


2006/11/21 


2454060 


55.0 




3G 


1996/11/22 


2450410 


47.0 




112 


2007/01/13 


2454114 


54.5 




37 


1997/01/06 


2450455 


47.6 


4 


112 


2007/01/14 


2454114 


55.0 




38 


1997/03/02 


2450510 


48.1 




113 


2007/03/10 


2454169 


54.9 




39 


1997/04/15 


2450554 


47.6 


4 


113 


2007/03/11 


2454171 


53.7 


6 


40 


1997/06/01 


2450601 


46.2 


4 


114 


2007/05/03 


2454224 


53.6 


6 


42 


1997/09/04 


2450696 


46.0 




118 


2007/12/03 


2454438 


52.9 




44 


1997/12/03 


2450786 


45.8 


3-4 


118 


2007/12/04 


2454438 


51.4 




45 


1998/01/16 


2450830 


44.8 


3-4 


119 


2008/01/20 


2454486 


45.6 




4G 


1998/03/05 


2450878 


44.9 


3-4 


120 


2008/03/03 


2454529 


46.0 


3-4 


47 


1998/04/15 


2450919 


46.0 




120 


2008/03/04 


2454530 


47.0 




48 


1998/06/01 


2450966 


49.7 


5 


121 


2008/04/21 


2454578 


49.3 




49 


1998/07/21 


2451016 


52.2 




121 


2008/04/23 


2454580 


49.2 


3-4 


50 


1998/09/21 


2451078 


52.8 




122 


2008/06/10 


2454628 


49.7 




51 


1998/11/06 


2451124 


50.6 




122 


2008/06/11 


2454629 


49.6 




52 


1998/12/28 


2451176 


49.2 




125 


2008/11/07 


2454778 


50.1 


3-4 


53 


1999/02/10 


2451220 


49.8 




125 


2008/11/08 


2454778 


49.8 




54 


1999/04/07 


2451276 


51.6 


>4 


126 


2008/12/29 


2454830 


49.6 




55 


1999/05/24 


2451323 


52.6 


5 


127 


2009/02/14 


2454877 


48.9 




57 


1999/09/11 


2451433 


51.0 




127 


2009/02/15 


2454878 


48.2 


4 


58 


1999/10/31 


2451483 


51.0 




128 


2009/04/04 


2454926 


51.2 




59 


1999/12/13 


2451526 


47.6 


4-5 


128 


2009/04/05 


2454926 


47.4 


4 


GO 


2000/02/06 


2451581 


46.0 




129 


2009/05/29 


2454981 


47.0 


5 


61 


2000/03/18 


2451622 


46.2 


5 


134 


2010/01/14 


2455211 


46.7 




G2 


2000/04/30 


2451665 


44.7 


5 


135 


2010/03/01 


2455257 


45.8 




G3 


2000/06/19 


2451715 


46.7 


>5 


135 


2010/03/01 


2455257 


46.6 




64 


2000/08/01 


2451758 


47.0 


>5 


136 


2010/04/14 


2455301 


47.7 


4 
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Table 1. (Continued.) 



E 


Date 


JD 


Ts (d) 


iV NO 


E 


Date 


JD 


Ts (d) 


^Vno 


65 


2000/09/21 


2451809 


46.5 




136 


2010/04/20 


2455307 


46.0 




66 


2000/11/04 


2451853 


47.9 




137 


2010/06/05 


2455353 


45.6 


3-4 


G7 


2000/12/21 


2451900 


47.3 




141 


2010/12/03 


2455534 


45.2 




68 


2001/02/11 


2451952 


48.4 


>5 


142 


2011/01/16 


2455578 


45.2 




69 


2001/03/27 


2451996 


48.4 




143 


2011/03/04 


2455624 


45.8 


4 


70 


2001/05/17 


2452047 


48.6 


6 


143 


2011/03/04 


2455625 


48.1 


4 


74 


2001/11/26 


2452240 


48.6 




143 


2011/03/05 


2455626 


49.5 


4 


75 


2002/01/17 


2452292 


48.0 


6 


144 


2011/04/22 


2455674 


49.5 


4 


76 


2002/03/01 


2452335 


45.0 




144 


2011/04/23 


2455675 


49.0 


4 


77 


2002/04/18 


2452383 


45.1 


>4 


144 


2011/04/24 


2455676 


48.9 


4 


81 


2002/10/10 


2452558 


46.1 




145 


2011/06/11 


2455724 


48.3 




82 


2002/12/02 


2452611 


46.6 




145 


2011/06/11 


2455724 


48.9 




83 


2003/01/23 


2452663 


50.0 




147 


2011/09/16 


2455820 


50.6 




84 


2003/03/10 


2452709 


49.2 


>5 


148 


2011/11/07 


2455872 


51.3 




85 


2003/04/29 


2452759 


48.3 


6-7 


149 


2012/01/02 


2455929 


53.0 




85 


2003/04/30 


2452760 


48.2 




149 


2012/01/02 


2455929 


53.9 




86 


2003/06/16 


2452807 


47.9 




150 


2012/02/20 


2455978 


53.4 




88 


2003/09/20 


2452903 


45.9 




150 


2012/02/27 


2455984 


53.8 




88 


2003/09/21 


2452903 


45.1 




152 


2012/06/11 


2456089 


54.4 




89 


2003/10/27 


2452940 


43.5 




152 


2012/06/12 


2456090 


52.8 





